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ABSTRACT 

Context. Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel 

space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation 

conditions, and far-IR photometry to constrain the dust distribution. 

Aims. We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. 

Mettiods. We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at A ~60-190 /im using the Photodetector 

Array Camera and Spectrometer instrument on the Herschel space observatory. 

Results. We detect [01] 63 fim emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, 

including [Cll] 158 and [Ol] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 yum. Extensive models are 

presented in associated papers. 

Key words. Stars: pre-main sequence; (Stars:) planetary systems: protoplanetary discs; Infrared: stars 
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1. Introduction 

Circumstellar discs of gas and dust are the environments in 
which planets form; thus, studies of the evolution of these discs 
will inform our understanding of the formation of gas giants. The 
' Herschel space observatory (Pilbratt et al. 1201 Oi l, with its facil- 
' ities for sensitive spectrometry and photometry, allows for sig- 
nificantly deeper searches for circumstellar gas than were previ- 
ously possible, providing tighter constraints on the timescale of 
. giant planet formation and the mechanisms of circumstellar disc 
■ evolution. 

The Gas in Protoplanetary Systems Herschel open time key 
project (GASPS, PI. W.R.F. Dent, see Dent et al., in prep.) aims 
' to characterize the evolution of gas and dust in circumstellar 
, discs. Using the Photodetector Array Camera and Spectrometer 
(PACS, Poglitsch et al. 2010) integral field unit spectrometer and 
bolometer array, we are carrying out a 400 hour survey of 240 
young stars for gas line emission and sensitive far-IR photome- 
try. We focus on three primary emission lines, [OI] 63.184 fim, 
[OI] 145.525 yum, and [CII] 157.741 jjm, as well as several emis- 
sion lines of H2O and CO and broadband photometry at 70 and 
160 fim to measure dust emission. Our survey covers young 
stars at ages from 1-30 Myr, in the Taurus, Upper Scorpius, rj 
Chamaeleontis, TWA, Beta Pictoris, and Tucana Horologium as- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



sociations, as well as HAeBe stars at a variety of ages 0.1- 30 
Myr, and with masses ranging from -0.5^ Mq. We chose sys- 
tems known from near-IR and millimeter photometry and accre- 
tion indicators to have discs in a variety of states, from those 
still accreting to debris discs, as well as a sample of stars show- 
ing no signs of a disc. We focus on optically-visible stars with 
low extinction and small envelope masses, i.e. SED Class II-III. 

At a basic level, observations of [CII] and [OI] will provide 
a determination of the presence of gas. Comparison to photo- 
dissociation models (e.g. Kaufman et al. 1999) allows for the 
use of [CII] and [OI] emission to place constraints on the density 
and temperature structure of circumstellar discs. Such observa- 
tions will help illuminate disparities between disc mass estimates 
based on sub-mm continuum measurements of dust emission and 
gas tracers such as CO. However, atomic C and O have emis- 
sion lines in the far-IR that are largely unobservable from the 
ground. While the Infrared Space Observatory (ISO) allowed 
for studies of atomic emission from a few bright sources (e.g. 
Lorenzetti et al. 1999), the Herschel PACS instrument has 2 or- 
ders of magnitude better point-source sensitivity than the ISO 
Long Wavelength Spectrograph, opening up the study of circum- 
stellar disc gas to hundreds more sources. 

To aid with the interpretation of our PACS line and con- 
tinuum data, we have interfaced tools for calculating dust and 
line radiative transfer (MCFOST, Pinte et al. 120061 ) and gas 
thermo-chemistry (ProDiMo, Woitke et al. 120091 ). We have used 
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these codes to produce a grid of 300,000 disc models, the Disc 
Evolution with Neat Theory (DENT, Woitke et al. l20TUb . 

In this paper, we present the first results of science demon- 
stration phase observations of a Herbig Ae star (HD 169142), 
two T-Tauri stars (RECX 15 in ?; Cha and TW Hya), and a debris 
disc system (HD 181327). All have ages in the range 6-12Myr, 
and bear discs that have undergone some evolution. Near/mid- 
IR photometry and mm photometry have indicated disc clearing 
or grain growth in the inner discs of HD 169142 and TW Hya 
(e.g. Grady et al. 120071 Hughes et al.'2007'). Mass accretion rate 
estimates range from < 7 x 10"'° for HD169142 (Grady et al. 
12007 ) to 5 X 10-1° - 2 X 10-'' Mo/yr for TW Hya (Muzerolle et 
al-EOOOi Her czeg et al. l2004i : and -lO'* Mo/yr for RECX 15 
(Lawsonetal. llfel i. 

In Sect. 121 we describe our observations and data reduction, 
and present our results. We present a preliminary discussion of 
the implications for the states of these discs in Sect. [3] 



2. Observations and results 

For our SDP targets, we carried out both spectroscopy and pho- 
tometry with the PACS instrument on Herschel. For each tar- 
get, we carried out a 1669s PacsLineSpec observation, a 5150s 
PacsRangeSpec observation, and either a 133s PacsPhoto scan 
map (RECX 15 and TW Hya) or 159s PacsPhoto point source 
observation (HD 169142 and HD 181327). We carried out an 
additional 220s scan map observation of RECX 15, to test obser- 
vation modes. We carried out the reduction of our observations 
with the facility reduction tool, HIPE (Ott, S. '2010), using the 
scripts disseminated at the Herschel data reduction workshop in 
January, 20100 



2.1. Spectroscopy 

Spectroscopic observations with the PACS IFU were carried out 
in chop-nod mode in order to remove the telescope emission. 
The PACS IFU is a 5x5 filled array, with 9.4"x9.4" spaxels. 
Observations were carried out such that the central IFU spaxel 
was centered on the target in both nods. Our SDP observa- 
tions were carried out with a 2" dither Each PacsLineSpec ob- 
servation includes simultaneous observations from 62.93-64.43 
and 180.76-190.29 fim. Our PacsRangeSpec observations in- 
cluded, in series, simultaneous observations from 71.81-73.28 
and 143.59-146.53 /urn, 78.37-79.76 and 156.70-159.47 pm, 
and 89.28-90.48 amd 178.51-180.96 |/m. Typical continuum 
rms for all regions are ~0.3 Jy, with Icr line flux uncertainties 
of 2-5x10"'** W/m^. Spectral resolution ranges from R~1100- 
3400. Full details of the observations are given in Dent et al. (in 
prep). 

For our point source flux measurements, we extracted the 
A and B nods for the central IFU spaxel. For each nod, we then 
Nyquist binned in wavelength with non-overlapping bins half the 
width of the instrumental resolution. The mean of the two nods 
restores the flux of the objects. We apply wavelength dependent 
flux and aperture corrections released by the PACS development 
tearrQ and assume a flux calibration uncertainty of 40%. 

To extract line fluxes and measure continuum flux from our 
reduced, nod-combined, flux corrected spectra, we carried out 
an inverse-error weighted fit of a gaussian and first-order poly- 
nomial to the spectra. Using the polynomial fit to each spectral 
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' https://nhscsci.ipac.caltech.edu/sc/index.php/Workshops/ 
^ PacsSpectroscopyPerformanceAndCalibration_v 1 _1 _0.pdf 



Fig. 1. Baseline subtracted, aperture and flux corrected [OI] 
63 fim spectra. Error bars indicate Icr statistical errors; gaussian 
line fits are shown in red. 



region, we estimated the continuum emission at the rest wave- 
length for each line. To estimate the error on the continuum, we 
calculated the error of the mean for the residual of the poly- 
nomial fit in a region from 2 to 10 instrumental FWHM from 
the line rest wavelength. For emission lines, we report the flux 
of detected lines as the integrated flux of the gaussian line fit. 
We calculate Icr fluxes as the integral of a gaussian with height 
equal to the continuum RMS, and width equal to the instrumen- 
tal FWHM; we report 3cr upper limits for non-detections. 

The [OI] 63 fim emission line is clearly detected towards HD 
169142, RECX 15, and TW Hya from our SDP sample; with 
FWHM -90 km/s, HD 169142 and TW Hya are spectrally unre- 
solved. The emission from RECX 15, on the other hand, appears 
much broader, with a FWHM ~150 km/s; this is a marginal re- 
sult and could be due to off-center or low S/N observations in 
the two nods (further discussion in Woitke et al., in prep.). We 
show our baseline subtracted [OI] 63 yum spectra in Figure [T] 
Examination of spectra from the off-center data cubes show no 
indications of extended or field emission. In Table [1] we show 
the [01] 63, [OI] 145, and [CII] 157 fluxes and 3cr upper Hmits, 
as well as the measured continuum fluxes and 3cr upper limits. 
At this stage, no other lines were detected with high significance. 



2.2. Photometry 

Simultaneous 70 and 160 fim continuum observations were ob- 
tained with the PACS photometer. The blue camera consists of 
a 64x32 array of 3.2" pixels, while the red camera consists of 
a 32x16 array of 6.4" pixels. Point source mode consists of 50" 
chod-nod observations for sky subtraction and a 3 point dither 
(8.5") for higher spatial sampling. Scan maps were carried out 
with a 20"/s scan rate, with a 10 Hz integration cycle. 

All sources appear spatially unresolved, with FWHM ~5" 
and 1 1" at 70 and 160 yum, respectively. We carried out photom- 
etry on images at the instrumental pixel size using a 21" aper- 
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Table 1. Stellar properties and observational results 



Object 



SpT 



Age 
Myr 



Distance 
pc 



[OI] 63° 
10-'* W/m- 



cont. 

Jy 



[01] 145° 
10""* W/m^ 



cont. 

Jy 



[CII] 157° 
10""* W/m- 



cont. 

Jy 



70 ^m' 

Jy 



160 Mi^f 

Jy 



HD 169142 


A5(" 


6+6(2) 


145") 


71.7±3.8 


18.89 


< 10.4 


13.76 


<6.4 


14.49 


HD 181327 


F5(4) 


12+8(5) 


50.6'^' 


<9.3 


1.75 


<8.0 


0.80 


<8.5 


0.77 


TWHya 


K7e<" 


1Q+10(8) 


56*" 


36.5±3.6 


2.99 


<5.5 


7.00 


<6.6 


8.79 


RECX 15 


M2e("" 


9+2(11) 


97(12) 


31.5±3.1 


0.14 


<4.6 


0.10 


<7.4 


<0.09 



27.35+0.03 17.39±0.05 

1.90+0.03 O.87±0.O3 

3.90+0.02 7.38±0.04 

0.18+0.02? 0.07±0.0is 



<"' Calibration uncertainty 40% <*> Itr ~0.03-0.05 Jy <^' Icr ~0.02-0.04 Jy <''> \cr ~0.03-0.05 Jy ''> Calibration uncertainty 5% <" Calibration 

uncertainty 10% *^' Results from two stacked scan maps. 



References. (1) D unkin e t al., [19971 (2) Gra dy et al ., [20071 (3) van Boekel, et al., [20051 (4) Nordstrom et al. [20041 : (5) Zuckerman & Song[2004[ : 
(6) Ferryman et al.TgW; 7) de la Reza et al. [19891 (8) Barrado Y Navascues [2056l (9) Wichmann et al. [19981 (10, 11) Lawson et al., l20D2l (12) 
Mamajek et al.. .2000. 



ture for both blue and red camera observations. The PACS team 
has reported flux corrections and aperture correction^ result- 
ing in a flux calibration uncertainty of 5% in the blue camera, 
and 10% in the red camera. All four of our SDP targets are de- 
tected in both bands; we report these continuum measurements 
without color coiTection in Table [T] None of the sources show 
signs of extended emission. Comparison with publicly available 
IRAS and Spitzer photometry shows a general agreement within 
the systematic uncertainties. Our model fitting efforts of lines 
and continuum from HD 169142 and TW Hy a usin g MCFOST 
and ProDiMo, are presented in Meeus et al. (120101 and Thi et 
al., (|20T0 ): RECX 15 and HD 181327 are addressed in further 
detail in Woitke et al. (in prep) and Menard et al. (in prep). 

We find one exception to a good agreement between contin- 
uum measurements from spectroscopy and photometry - blue 
camera observations of HD 169142 show differences up to 50%. 
The cause of the difference is currently unknown, and may be 
due to the large systematic flux uncertainties. We will revisit this 
issue with future revisions of the instrumental calibration and 
reduction techniques. As the 70 //m flux from the photometric 
observation is more consistent with photometry in the literature, 
we consider it the more reliable of the two flux values, and note 
the [OI] fine flux for HD 169142 may be higher by up to 50%. 



3. Discussion 

3.1. Discs or ouflows? 

A primary driver of the GASPS project is to use far-infrared 
emission lines to characterize the gas content of circumstellar 
discs. However, many young star with accreting discs are also 
known to drive outflows. Their shocks produce strong line emis- 
sion in this band (Liseau et al. 2006). Following HoUenbach 
(I1985I I. we can bound the mass outflow rate from the [OI] 
63 luminosity, Moutflow(M0/yr) ^ 10"''L[oi]63(Lo). As some 
shock emission may emerge through other lines, for the ob- 
served fluxes, this translates into lower limits of 46, 3.5, and 
9.1 X 10"' Mo/yr for HD 169142, TW Hya, and RECX 15 re- 
spectively. Comparison to observed accretion rates (Sect.[T]) rule 
out an outflow shock origin for the emission from HD 169142, 
and requires the outflow rate to be a very high fraction (> 0.2) 
of the accretion rate for TW Hya and RECX 15. This is much 
greater than the ratio found in T Tauri stars by Hartigan et al. 
([T995l l. 

Further, there is no direct evidence for outflows in these 
sources. We do not find [OI] in other IFU spaxels in excess of 



that expected in the wings of the point-spread function, which 
places an upper bound to the extent of the emission to less than 
680, 260, and 460 AU for HD169142, TW Hya, and RECX 15 
respectively. With the previously noted exception of RECX 15, 
the spectra are not resolved, precluding jets with line-of-sight ve- 
locities greater than ~45 km/s. As noted in Sect. 12.1! the broad 
profile in RECX 15 is suspected to be due to instrumental effects. 
This also agrees with the naiTow spectral profiles of OI 6300 A 
an d CO rotational lines in HD 169 142 and TW Hydra (Acke et 
al. [20051 Raman et al. .2006. Qi et al. .2004) . 



While J and C shock models predict 



[0/]63 



and 



[0/]63 



line 



PACS.ScanMap_ReleaseNote.23Feb2010.pdf 



[0/ 1145 "■■" [C//J158 

ratios >10 (e.g. HoUenbach & McKee [T989 '). the analysis of 500 
archival ISO observations by Liseau et al. (12006) showed typical 
ratios of ~ 1-20 and 0.1-5 for background subtracted sources (lo- 
cus of detections of Liseau et al., from their Fig. 2, are shown in 
Fig.|2] along with models from the DENT grid as discussed be- 
low). The high ^j^ ratios of our targets (>10 for HD 169142, 
and >5 for TW Hya and RECX 15) only marginally overlap the 
outflow locus - less than 10% of the Liseau et al. sources have 
line ratios large enough to be consistent with our detections. 

We compare our results with a subset of the DENT grid of 
models, selected to match the observed range of [OI] 63 fluxes 
(scaled to 140 pc). This includes all models with [OI] 63 line 
strengths from (4-90)xl0""^ W/m^. We show the locus of these 
models in Fig.|2l which include models with disc masses ranging 
from 10"^ to 10"' Mq. Location within this plot can generally 
constrain disc gas mass to within an order of magnitude; further 
constraints using detected line strengths, stellar properties and 
data from other sources such as optical and mm continuum and 
line measurements, can refine this estimate and allow for deter- 
mination of the properties of the originating structure (see Kamp 
et al.. l2010l and Woitke et al., 2010). Predicted line flux ratios of 
the best fit models for HD 169142 (Meeus et al. 12010) . TW 
Hya (Thi et al.,2010), and for RECX 15 (Woitke et al., in prep) 
indicate that with somewhat longer observations, TW Hya and 
RECX 15 may be detectable in [CII] 158, whereas HD 169142 
may be detectable in [OI] 145. 

Given the lack of high velocity spectral observations, poor 
overlap with observed line ratios of jet sources, and consistency 
of line ratios with disc models, we cautiously interpret our re- 
sults in the context of disc emission. 



3.2. [CII] non-detections 

The lack of detections of the [CII] 158 line from any of our SDP 
targets is surprising. Kamp et al. ( 12010! ) predict [CII] 158 emis- 
sion of -10"'^ W/m^ for HAe discs with radii -200 AU, with 
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Fig. 2. Line flux ratio diagram comparing models with observa- 
tions. Dots represent GASPS observations with Hnes indicating 
3cr limits, and arrows indicating lower limits. Squares show pre- 
dicted line flux ratios based on the detailed models (HD 169142: 
Meeus et al.. l20T0l TW Hya: Thi et al., 2010; RECX 15: Woitke 
et al., in prep.). Outflow sources are shown in red (Liseau et al. 
120061 lines indicate contour containing 90% of points, hash indi- 
cates 68% of points.). Gray regions indicate 99, 90, and 68% of 
DENT models with [OI] 63 flux consistent with our observations 
(4-90x10-1** W/m^ at 140 pc). 



emission varying with disc surface area. HD 169142 is known 
from mm imaging to have a disc of radius >200 AU (Raman et 
al. [2006). Other factors must be at play, such as flat discs, lack 
of UV flux, or shielding by a large amount of small dust grains 
at large radii. 



3.3. A new regime of disc studies 

The ratio of [OI] 63 line strength to continuum varies sub- 
stantially, from 0.3-1 for HD 169142 and TWHya to 17 for 
RECX 15. Both line and continuum are likely optically thick and 
the stellar luminosities are very diff'erent for the three sources 
but one possible explanation for the strikingly high line-to- 
continuum ratio in RECX 15 is an enhanced gas-to-dust ratio 
(see Fig. 3 of Pinte et al.. l20T0t . 

Further, whereas HD 169142 and TWHya have strong mm 
continuum fluxes indicating large dust masses, RECX 15 is un- 
detected at 850 /^m (Phillips, in prep.). This is similar to the 49 
Ceti disc which is detected at 1.3 mm in CO 2-1 but not in the 
continuum (Hughes et al. "2008). These comparisons suggest a 
high gas-to-dust ratio; however, further modeling is necessary 
(see Woitke et al., in prep.). The GASPS survey will reveal how 
common are discs with gas emission yet lacking detectable dust. 

Prior studies of circumstellar disc gas have focused on near- 
infrared emission lines from the hot inner disc (e.g. Brittain et 
al. 12007) or mm lines from the cool outer disc (e.g. Dent et al. 
I2005I I. Far-infrared line emission from intermediate radii of gas 
rich discs has not yet been explored in similar detail. Thanks 
to Herschel's unprecedented sensitivity, GASPS will examine 
the gas content of circumstellar discs in previously unaccessible 
regimes. 
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